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ABSTRACT 

A calculation of the total cross section for top quark production in hadron- 
hadron collisions is presented based on an all-orders perturbative resummation of 
initial-state gluon radiative contributions to the basic quantum chromodynamics 
subprocesses. Principal-value resummation is used to evaluate all relevant large 
threshold contributions. In this method there are no arbitrary infrared cutoffs, 
and the perturbative regime of applicability is well defined, two attributes that 
significantly reduce the estimated uncertainty of the results. For pp collisions 
at center-of-mass energy = 1.8 TeV and a top mass of 175 GeV, we obtain 
aitt) — 5.52lQ4gp6, in agreement with experiment. Predicted cross sections are 
provided as a function of top mass in pp collisions at ^/s = 2.0 TeV and in pp 
collisions at CERN LHC energies. 

1. Introduction 



The quest for the top quark t reached fruition in the past year with the observa- 
tion of ti pair production in proton-antiproton colhsions at the Fermilab TevatronB. 
A deserving question is the quantitative rehabihty of theoretical computations of the 
total cross section, as a function of top mass, based on the main production mecha- 
nisms in perturbative quantum chromodynamics (pQCD). In this paper, we discuss 
the motivation for incorporating the effects of initial-state gluon radiative corrections, 
and we present our all orders resummation of these contributions.! 

At lowest order in perturbation theory, two QCD partonic subprocesses contribute 
to p + p — > t + i + X. They are quark-antiquark annihilation: q + q ^ t + i and 
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Fig. 1. The parton-parton cross sections as a function of rj in the AIS scheme at m = 175 GeV for 
the subprocesses (a) qq — > tiX and (b) gg — > tiX. Plotted are the lowest order Born cross section 
(dotted line) and the next-to-leading order cross section (solid line). The QCD scale /i = m. 

gluon-gluon fusion: g + g t + i. Short-distance partonic cross sections based on 
these lowest order C(a^) subprocesses and on the next-to-leading 0{a^) subprocesses 
have been investigated thoroughly.!'! A full 0{a'^) calculation, for n > 4, does not 
exist. The physical top quark cross section is obtained from a convolution of the 
perturbative short-distance subprocess cross sections with parton distributions that 
specify the probability densities of the quarks, antiquarks, and gluons of the incident 
p and p. Our work addresses improvements in the reliability of calculations of the 
subprocess cross sections. 

The motivation for this work begins with the observation that the size of the 
0{al) terms in the qq and gg partonic cross sections are much larger than their 
0{a1) counterparts in some kinematic regions, notably in the near-threshold region 
of small rj. The variable 1] = s/Am? — 1, where s is the square of the energy of parton- 
parton subprocess and m denotes the mass of the top quark. Variable rj measures the 
"distance" above the partonic production threshold. 

As shown in Fig. 1, for a top mass of 175 GeV, in both the qq and the gg channels 
the size of the C(a^) term exceeds that of the 0{al) term for rj ~ 0.1, and the 
ratio grows as rj decreases.! Therefore, the important notion underlying perturbation 
theory, that successive terms in the perturbation series should be smaller, is not valid 
at small r/, i.e., in the region near production threshold. This region of phase space 
is important for top quark production at the Tevatron. Owing to the large mass of 
the top quark, relative to the pp center of mass energy ^/s, the near-threshold region 
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contributes significantly when the convolution integral, mentioned above, is done over 
the full of range of r]. Confidence in the results of a perturbative calculation of the 
overall ti cross section requires an appropriate understanding of the origin of the large 
next-to-leading order enhancement of the partonic cross sections near threshold. (In 
the gg channel, the ratio of the 0{a^) and C(a^) terms exceeds unity for large r] also. 
The gg channel and the large rj region are important for bottom quark production at 
the Tevatroni, but not for top quark production.) 

We treat only ti pair production. Mechanisms for single top productioni are not 
considered here. At m = 175 GeV and y/s = 1.8 TeV, they contribute a cross section 
about 20% that of the pair production mechanisms. 

2. Gluon Radiation and Resummation 

The origin of the large threshold enhancements in the subprocess cross sections 
may be traced to initial-state gluon radiation.i After the cancellation of soft singular- 
ities between the contributions from real gluon emission and virtual gluon exchange 
terms, and proper factorization of coUinear divergences, there remain terms at 0{al) 
that are proportional to in{l — z). The variable z = 1 — 2kg.pt/'m'^ where pt and kg 
are the four-vector momenta of the produced top quark and the gluon radiated into 
the final state in the 2-^3 process. The limit z —>■ 1 corresponds to zero momentum 
carried by the gluon. 

The partonic cross section may be expressed generally as 



crij{ri,m ) = I dz 
In the near threshold region. 



Hijiz, a) ~ 2aCijin\l - z) + 



1 + Hij^z, a] 



Kjiv^m^.z). (1) 



2Cl/n\l -z)- \CiJh2^n\l - z) 



(2) 



We work in the MS factorization scheme in which the g, q and g densities and the 
next-to-leading order partonic cross sections are defined unambiguously. In Eq. (|I]), 
the lower hmit of integration Zmin = 1 — 4(1 + rf) + 4a/1 + 77, and ij E {qq,gg} 
denotes the initial parton channel. We set a = as{m)/TT. Symbol a[j{ri^im? , z) = 

d{d'if {r],m'^, z))/dz, where a^f is the lowest order partonic cross section expressed 
in terms of inelastic kinematic variable^ to account for the emitted radiation. The 
integration in Eq. (|l|) is over the phase space of the radiated gluons, parametrized 
through the dimensionless variable z. In Eq. (^, Cij is the color factor for the ij 
production channel. 

Equation (0) approximates the near-threshold behavior of the partonic cross sec- 
tion. It manifests the logarithmic behavior in{l — z) mentioned above. Explicit 
calculationJi of the complete 0{a^) cross section provide the 2aCij£n^{l — z) term. 
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The terms proportional to o? are appropriated from 0(0^) computations of massive 
lepton-pair pro duct ion based on the assumption of universahty of leading loga- 
rithmic contributions. As in other hard-scattering processes, where large logarithmic 
contributions are present near threshold, the goal of gluon resummation in tt pro- 
duction is to sum the series in a"£n^"(l — z) to all orders. Resummation, studied 
extensively for massive lepton-pair production, i'ilisl is important both for theoret- 
ical understanding of the perturbative process and for stability of the quantitative 
predictions. 

In resummation procedures, the large logarithmic contributions are exponentiated 
into a function of the QCD running coupling evaluated at a variable momentum scale 
that is a measure of the radiated gluon momentum. A straightforward method of 
resummation for tt production was published a few years ago.i In this approach, 
the partonic cross section of Eq. (|lD is replaced in the MS scheme by the resummed 
expression 



where 



E^f^iz, m^) oc Ci,a{{l - zf/V)in\l - z). (4) 



It is easy to verify that the form of Eqs. (|T|) and (0) is reproduced if e^'j is expanded 
in a power series in a(m^). A limitation of this method is that an infrared singularity 
is encountered in the soft-gluon limit z ^ \: owing to the logarithmic behavior of 
a(g^), Q;(g^) oc ln~^{q^ / !s?q(jj^), a((l — zY^'^m?) oo as z ^ 1. This divergence 
of the integrand at the upper limit of integration necessitates introduction of the 
undetermined infrared cutoff cutoff /Xq in Eq. (0), Aqcd < /Wo < that serves to 
prevent the integration over z from reaching the Landau pole of the QCD running 
coupling constant. The cutoff has a related effect of eliminating a portion of the 
integration over the partonic subenergy when the convolution with parton densities 
is done to obtain the physical cross section. In terms of 1], the convolution is restricted 
to 1] > r]o = (/io/'^)^/2. The presence of an extra scale spoils the renormalization 
group properties of the overall expression. Moreover, dependence of the resummed 
cross section on this undetermined cutoff is important numerically.! 

Laenen et al furnish their final predictions in the DIS factorization scheme in 
which Eqs. (|) and (^) are modified slightly. They obtain 

au{m = 175 GeV) = 4.95+[!:^ pb , (5) 

based on the assumed values fio = 0.1m for the qq channel and fio = 0.25m for the 
gg channel. Owing to sensitivity to /io, it is difficult to assess the significance of the 
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estimated uncertainties. 



3. Principal Value Resummation 

The principal-value method of resummation (PVR)i has an important technical 
advantage in that it does not depend on arbitrary infrared cutoffs, as all Landau-pole 
singularities are by-passed by a Cauchy principal-value prescription. Because extra 
undetermined scales are absent, the method also permits an evaluation of the pertur- 
bative regime of applicability of the method, i.e., the region of the gluon radiation 
phase space where perturbation theory should be valid. The method has been tested 
successfully in massive lepton-pair production.0 

To illustrate how infrared cutoffs are avoided in the PVR method, it is useful to 
express in moment space the exponent that resums the in{l — z) terms: 

i?(n,m^) = — y dx"^—^ j -^g a (^Xrn^^ . (6) 

(l-x)2 

The function g{a) is calculable perturbatively, but, again, the behavior of a(Am^) 
leads to divergence of the integral when Am^ — > Ag^^,. To tame the divergence, a 
cutoff can be introduced in the integral over x or directly in momentum space, in 
the fashion of Laenen et al .0 In the principal-value redefinition of resummation, the 
singularity is avoided by replacement of the integral over the real axis x in Eq. (|^) by 
an integral along a contour P in the complex plane: 

E^>,m^)^-/rfC^^^ / y^7[«(Am^)]. (7) 

P ^ (1-C)2 

The function E^^ {n,m'^) is finite since the Landau pole singularity is by-passed. In 
Eq. (^, all large soft-gluon threshold contributions are included through the two- loop 
running of a. 

Equations (|^) and (|^ have identical perturbative content, but, when expanded 
in power series in a(m^) and in Aqcd/^, they manifest differences in their inverse 
power (high-twist) terms. Since the inverse power content is not a prediction of 
perturbative QCD, neither expression is a priori preferable, except for the attractive 
finiteness of Eq. (^. In our application of principal- value resummation to top quark 
production, we choose to use the result only in the region of phase space in which 
the perturbative content dominates. Thus, the high-twist content of Eq. (^), and the 
difference between the high- twist components of Eqs. (|^) and @, are not matters of 
phenomenological significance. 

After inversion of the Mellin transform, the resummed partonic cross sections 
according to PVR, including all large threshold corrections, can be written in the 
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form of Eq. (|1|) , but with Eq. (0) replaced by 

/•ln(— ^) °° 

'H^y[z,a)= rfxe^-(^'")5]Q,(x,«) . (8) 

The leading large threshold corrections are contained in the exponent Eij{x, a), a cal- 
culable polynomial in x. The functions {Qj{x, a)} arise from the analytical inversion 
of the Mellin transform from moment space to the physically relevant momentum 
space expressed in Eq. (|). These functions are produced by the resummation and 
are expressed in terms of successive derivatives of E: Pk{x, a) = d^E{x, a)/k\d^x. 

The functional form of Eij for tt production is identical to that for // production, 
except for the identification of the two separate channels, denoted by the subscript 
ij. However, only the leading threshold corrections are universal. Final-state gluon 
radiation as well as initial-state/final-state interference effects produce sub-leading 
logarithmic contributions that differ for processes with different final states. Among 
all {Qj} in Eq. (|^), only the very leading one is universal. This is the linear term 
in Pi, which turns out to be Pi itself. Since we intend to resum only the universal 
leading logarithms, we retain only Pi. Hence, Eq. (|]) can be integrated explicitly, 
and the resummed version of Eq. is 

ni^) = r™"' d;2e^-(''^(T^)'")<j;.(r/, z). (9) 

The upper limit of integration in Eq. is set by the boundary between the per- 
turbative and high-twist regimes. To characterize a region in moment space as high- 
twist, one must convert to momentum space through inversion of the Mellin trans- 
form, Eq. (^. Specification of the boundary is realized by the constraint that all 
{Qj}i J ^ 1 be small compared to Qq. This constraint can be shown to correspond 
to 

Pi(^ln(^^),a) <1. (10) 

As remarked above, we accept only the perturbative content of principal- value resum- 
mation, and our cross section is evaluated accordingly. Specifically, we use Eq. (^) 
with the upper limit of integration, Zmaxi calculated from Eq. (|TUp. The upshot is an 
effective threshold cutoff on the integral over the scaled subenergy variable ?7, reminis- 
cent of that introduced by Laenen et al , but one that is calculable, not arbitrary. In 
our case, the cutoff restricts the region of applicability of resummation to the part of 
phase space in which the perturbative content of Eq. (H) is the dominant content. For 
the top mass m = 175 GeV, we determine that the perturbative regime is restricted 
io 7] > 0.007 for the qq channel and rj > 0.05 for the gg channel. The difference 
reflects the larger color factor in the gg case. (One could attempt to apply Eq. (||) 
all the way to z^ax = 1, i-e., to ?7 = 0, beyond the perturbative regime of Eq. (p!OD, 
but one would then be using a model for non-perturbative effects, the one suggested 
by PVR, far beyond the knowledge justified by perturbation theory.) 
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Fig. 2. The parton-parton cross sections as a function of rj in the MS scheme at m = 175 GeV for 
the subprocesses (a) qq ttX and (b) gg tiX. Plotted are the lowest order Born cross section 
(dotted Hne), the next-to- leading order cross section (dashed hne), and the resummed cross section 
(solid line). The QCD scale fj, = m. 

Our final result does not rely critically on the PVR method to by-pass infrared 
renormalons and associated problems, precisely because we restrict application to 
the perturbative regime. In this regard, the presence of arbitrary infrared cutoffs 
in other resummation methods is superfluous, as all necessary information about 
infrared sensitivity (i.e., the perturbative regime) can be obtained by examining the 
perturbative asymptotic properties of the resummation functions. 

The resummation procedure includes only the leading threshold in'^{l — z) piece 
of the full calculation. To restore the full content of the complete next-to- 

leading order calculation, af^^^^ , we define our final resummed cross sections for each 
production channel through the improved prediction 

af^^\r^, m') = ^^-^(ry, m') + af^^'^ir^, m') - ^^(r^, m')\^^ . (11) 

The last term in Eq. ( [TTl) is the part of the next-to-leading order partonic cross sec- 
tion included in the resummation. In Fig. 2 we present the resummed partonic cross 
sections in the qq and gg channels at m = 175 GeV. We also show the lowest order 
and next-to-leading order counterparts. The three curves differ substantially in the 
partonic threshold region rj < 1, with the final resummed curve exceeding the other 
two. Below rj ^ 0.007 in the qq channel and rj ~ 0.05 in the gg channel, our resummed 
cross sections become identical to the next-to-leading order cross sections, a conse- 
quence of our decision to restrict the resummation to the perturbative domain. Above 
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P t f+X, v's= 1 ,5 TeV 



140 160 180 200 220 240 260 

m (GeV) 

Fig. 3. Physical cross section foiiJ^p — > {tt)X at y/s = 1.8 TeV as a function of top mass. Data from 
the CDF and DO coUaborationsEl are plotted. Shown are calculations for three choices of the scale 
^/m — 0.5 (dashed), 1 (solid), and 2 (dotted). 



77 ~ 1, our resummed cross sections are essentially identical to the next-to-leading 
order cross sections, as is to be expected since the near-threshold enhancements that 
concern us in this paper are not relevant at large 77. 

4. Calculations at ^ = 1.8 TeV and ^/s = 2 TeV 

In the remainder of this report, we present our results for the physical inclusive 
total cross section for ti production for a top- mass range m G {150, 250} GeV, includ- 
ing a discussion of the remaining theoretical uncertainties. We obtain the physical 
cross section by convoluting Eq. ([Tl|) with CTEQ3M parton densities^ and adding 
the results from both channels. In Fig. 3 we show the top-mass dependence of the 
physical cross section for pp — ^ (tt)X at =1.8 TeV. 

As illustrated in Fig. 4, the behavior of the physical cross section as a function 
of the renormalization/factorization scale fi is mild in the range /i/m G {0.5,2}. 
(The next-to- leading order result in Fig. 4 differs somewhat from that shown in 
our published paper.B Owing to a compiler error, the next-to-leading order curve in 
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Fig. 4. Plot showing the calculated dependence of the final resummed cross section on (/i/m) for 
tt production at m = 175 GeV and ^/s = 1.8 TeV. Shown also is the next-to- leading order result 
(dashed curve). 

the pubhshed figure is incorrect.) We consider the variation of the cross section in 
the range /x/m G {0.5,2} to be a reasonable measure of the theoretical perturbative 
uncertainty. Over this range, the band of variation of the strong coupling strength 
as is a generous ±10% . We determine 

a'\m = 175 GeV, = 1.8 TeV) = 5.52^°;°^ pb . (12) 

We define the central value (5.52 pb) to be that obtained with /x/m = 1. The upper 
and lower limits correspond to the maximum and minimum values of the cross section 
in the range fi/m G {0.5,2}. The cross section is insensitive to the choice of parton 
densities. Repeating the same analysis with the MRS (A') densities^, we obtain 

a"(m = 175 GeV, ^/s = 1.8 TeV) = 5.32^°;°^ pb . (13) 

The central values in Eqs. (|I2|) and (|T3|) are about 10% larger than that of Laenen 
et al , Eq. (|^), but within the quoted uncertainties. Our calculated cross sections fall 
within the current uncertainty bands of the CDF and DO experiments.0 

The bands of perturbative uncertainty quoted in Eqs. (|l^) and (^) are rela- 
tively narrow. On the other hand, we noted in discussing Eq. (|]) that 2;max < 1, 
meaning that there is a reasonable range of t] near threshold in which perturbative 
resummation does not apply. Perturbation theory is not justified in this region. Cor- 
respondingly, further strong interaction enhancements of the ti cross section may 
arise from physics in this region. We know of no reliable way to estimate the size of 
such non-perturbative effects and, therefore, cannot include such uncertainties in the 
estimates of the perturbative uncertainty of Eqs. ([T2|) and (p^). 
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p t t+X, Vs=2.0 TeV 
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m (GeV) 

Fig. 5. Physical cross section for pp — > {tt)X at ^/s = 2 TeV as a function of top mass. Shown are 
calculations for three choices of the scale ii/m = 0.5 (dashed), 1 (solid), and 2 (dotted). 

In Fig. 5 we show the top-mass dependence of the physical cross section for 
pp ~^ {tt)X at the shghtly larger energy i/i = 2 TeV. We predict 

a"(m = 175 GeV, y/s = 2 TeV) = 7.56;g;^^ pb . (14) 

At m — 175 GeV, the value of the cross section at a/s — 2 TeV is about 37% greater 
than that at = 1.8 TeV. 

5. LHC Predictions 

Turning to pp scattering at the energies of the Large Hadron CoUider (LHC) at 
CERN, we note a few significant differences from pp scattering at the energy of the 
Fermilab Tcvatron. The dominance of the qq production channel at the Tevatron is 
replaced by gg dominance at the LHC. Owing to the much larger value of ^/s, the 
near-threshold region in the subenergy variable is relatively less important, reducing 
the significance of initial-state soft gluon radiation. Lastly, physics in the region of 
large V§, where straightforward next-to-leading order QCD is also inadequate, may 
become significant for ti production at LHC energies. Using the approach described 
in this paper, focussed on the resummation of initial-state gluon radiation, we present 
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predictions in Fig. 6 for LHC energies of 10 and 14 TeV. We estimate 

cT**(m = 175 GeV, ^/s = U TeV) = 760 pb . (15) 




P t i+X, Vs=10 and 1 4 TeV 
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Fig. 6. Cross section for pp {tt)X at -^i =10 and 14 TeV as a function of top mass. Shown are 
the next-to- leading order (dashed) and resummed (sohd) calculations for the scale choice /u/m = 1. 

6. Summary and Discussion 

In summary, we present a calculation of the total cross section for tt pair pro- 
duction in perturbative QCD including resummation of initial-sate gluon radiation 
to all orders in as- Two advantages of the principal- value method of resummation 
are the well-defined perturbative domain of applicability and the absence of arbitrary 
infrared cutoffs. Both qq and gg production channels are included in the calculation. 
At ^/s = 1.8 TeV, our final resummed cross section is approximately 10% greater 
than the pure next-to-leading order result and in agreement with data. We provide 
predictions for the cross section as a function of top mass at a/s — 2, 10, and 14 TeV. 

We remark that the resummation method developed here and applied to tt pair 
production is relevant in other situations in which dynamics probes the near-threshold 
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region in the scaled subenergy variable. An example of current interest at ^/s =1.8 
TeV is the production of hadronic jets that carry large transverse momentum.lii 
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